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ABSTRACT 

A  theoretical  scheme  to  describe  the  infrared  contrast  field  left  by 
a  surface  ship  on  a  body  of  water  is  presented.  Starting  from  a  field 
velocity  distribution  behind  the  vessel,  the  fluid  mechanical  behavior  of 
the  disturbance  is  developed  as  it  interacts  with  the  vertical  thermal  pro¬ 
file  of  the  water  and  produces  a  temperature  pattern  on  the  sea  surface. 

As  a  result,  the  signal  sensed  by  a  remote  detection  system  is  obtained. 

A  detailed  sensitivity  study  of  the  computational  results  is  con¬ 
ducted  by  separately  varying  each  of  the  main  code  variables.  This 
method  yields  program  input  values  that  best  fit  the  measurements. 

Using  a  bell-shaped  velocity  distribution  as  an  input,  the  contrast 
radiant  intensity  of  a  natural  macrowake  generated  by  a  catamaran  is 
computed  to  be  -225  kW/sr.  This  compares  with  a  measured  value  of 
—298  kW/sr,  from  which  it  differs  by  about  24  % . 

The  variation  of  the  spatial  width  of  the  macrowake  on  the  water 
surface  is  computed,  and  found  to  fit  the  measurements.  An  infrared  gap 
is  sometimes  observed  at  the  onset  of  the  macrowake.  Actual  imagery  of 
a  catamaran  and  a  monohull  ship  shows  definite  evidence  of  gaps.  Gap 
extents  can  be  predicted  by  the  theory  to  better  than  5%  accuracy. 
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INTRODUCTION 

Whenever  a  surface  ship  moves  across  a  body  of  water,  it  leaves  behind  it  a 
characteristic  pattern  that  can  be  observed  for  many  kilometers.  Detection  is  achieved 
through  the  application  of  a  number  of  techniques,  including:  visible  photography,  acoustics, 
conventional  radar,  synthetic  aperture  radar,  microwave  radiometry  and  infrared  sensing, 
which  is  the  subject  of  the  present  report. 

Four  distinct  regions  characterize  the  macrowake.  Just  behind  the  ship,  it  begins  as 
foamy,  aerated  and  turbulent  water.  The  bow  wave  originated  by  ship  motion,  the  turbulent 
boundary  layer  adjacent  to  the  hull,  and  the  effect  of  the  propeller  all  produce  large  amounts 
of  foam. 

Following  the  initial  spreading  is  a  region  known  as  the  turbulent  core.  In  this  region, 
no  significant  production  of  white  water  occurs,  although  many  bubbles  are  still  to  be  found 
at  and  below  the  surface.  If  the  upper  layer  of  the  sea  is  thermally  stratified,  then  the  violent 
mixing  that  takes  place  in  the  turbulent  core  will  bring  to  the  surface  some  water  which  is  at 
a  different  temperature.  The  temperature  difference  between  the  fluid  inside  the  contrast  field 
and  the  surrounding  water  makes  the  macrowake  detectable  by  an  infrared  radiometer.  After 
the  disappearance  of  the  foam,  the  turbulent  core  remains  visually  distinct  from  the  surroun¬ 
ding  ocean  surface  because  of  its  slick  appearance.  This  second  region  ends  when  iheie  is  no 
longer  a  visual  distinction  between  the  ocean  surface  and  the  macrowake. 


Small  hubbies  still  remain  in  the  thermally -mixed  water  of  the  third,  region,  and  the 
temperature  remains  different  from  that  of  the  surrounding  water. 

In  the  final  stage,  bubbles  in  the  turbulent  core  have  disappeared,  but  the  thermally 
mixed  water  at  the  surface  has  not  returned  to  the  state  of  thermal  equilibrium  that  existed 
before  the  passage  of  the  ship.  When  it  does,  the  contrast  field  can  no  longer  be  detected. 

Some  pioneering  work  in  wake  analysis  can  be  traced  back  to  the  Second  World 
War.1  Unfortunately,  most  measurements  performed  since  then  lack  instrument  calibration 
as  well  as  any  sea-truth  information.  On  the  theoretical  side,  a  thorough  search  of  the 
literature  has  failed  to  uncover  any  previous  attempt  at  a  mathematical  understanding  of  the 
infrared  signal  generated  by  the  inacrowake  of  surface  ships. 

The  surface  manifestation  of  the  macrowake  starts  out  as  a  three-dimensional  flow 
pattern  within  the  water.  This  situation  is  conceptually  depicted  in  Fig.  1 .  When  real 
measurements  are  simulated  with  a  computer,  the  flow  can  be  made  to  interact  with  a  ver¬ 
tically  stratified  ocean.  The  in-band  contrast  radiant  intensity  of  the  macrowake  is  then 
calculated  at  various  distances.  The  close  range  value  is  compared  with  the  measured  data. 
The  spatial  width  of  the  macrowake  is  compared  to  the  corresponding  image  data.  In  addi¬ 
tion  to  presenting  the  details  of  an  actual  simulation,  this  paper  concludes  with  a  calcula¬ 
tion  of  the  thermal  gap  that  is  sometimes  observed  behind  surface  ships. 

METHOD 

The  main  features  of  interest  are  the  temperature  profile,  the  spatial  width  of  the 
macrowake  on  the  free  surface,  and  the  contrast  signal  generated  by  the  infrared 
macrowake  at  a  remote  sensor  system.  All  these  quantities  can  be  readily  compared  with 
accurate  measurements. 

The  principal  measurements  used  in  the  present  work  were  gathered  during  the  Phelps 
Bank  microwave  experiment,  which  was  conducted  to  investigate  the  anomalous  observa¬ 
tions  made  by  the  SEASAT  satellite.2,3  Mr.  P.  M.  Smith  of  the  Naval  Ocean  Research 
and  Development  Activity,  kindly  provided  a  number  of  digitized  infrared  images  of  the 
oceanographic  research  ship  Hayes  and  its  macrowake.  Mr.  Smith  also  sent  us  an  image  of 
the  macrowake  generated  behind  the  surveying  ship  Kane,  taken  a  few  years  earlier  in  the 
Sargasso  Sea.4 

The  details  of  the  analysis  can  be  followed  in  the  block  diagram  of  Fig.  2.  In  this 
analysis,  velocity  data  from  a  vertical  plane  behind  the  propeller  and  vertical  water  density 
information  are  read  by  the  computer.  A  fluid  mechanical  code  propagates  the  flow  in 
three  dimensions  within  the  water.  This  method  eventually  generates  a  temperature  map  of 
the  water.  The  pattern  on  the  free  surface  can  be  translated  by  the  detector  module  into  a 
contrast  signal.  The  inclusion  of  the  atmospheric  propagation  yields  the  infrared  signal  as 
it  is  received  by  a  remote  sensor. 

A  subsurface  wake  program  called  GORWAK5  was  used  for  the  fluid  mechanical 
code.  This  code  solves  the  inviscid  two-dimensional  equations  for  the  evolution  of  a  densi¬ 
ty  and  flow  disturbance  in  a  stratified  ocean  environment.  It  uses  a  finite-difference 
numerical  integration  of  the  Euler  equations  in  the  vorticity-stream-function  representation. 
The  GORWAK  description  includes  the  surface  and  bottom  of  the  ocean. 

Despite  obvious  advantages,  the  use  of  this  program  also  involves  possible  problems. 
One  may  question  the  suitability  of  using  a  subsurface  code  in  modeling  a  surface  contrast 


field  as  is  the  case  in  the  infrared.  Another  concern  is  that  an  underwater  wake  is  essen¬ 
tially  morncntumless,  while  a  surface  ship  generates  an  “ordinary”  wake,  which  is  one 
where  mixing  and  air  entrainment  occur. 

SIMULATION 

Consider  a  situation  in  which  the  sea  is  described  in  terms  of  a  number  of  layers. 
Displacing  a  vertical  column  of  the  fluid  infinitesimally  from  its  position  of  equilibrium 
gives  rise  to  a  natural  frequency  of  oscillation.  This  frequency  is  called  the  Brunt-Vaisala 
frequency6  and  is  defined  by  the  expression 


where 

N 

g 

z 

Co 

Q 


(If) 


0) 


Brunt-Vaisala  frequency  in  rad/s 
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Equation  1  can  now  be  rewritten  as 
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where  T  is  the  temperature  in  °C.  Since  the  thermal  expansivity  a  is  given  by 
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Density  values  for  ordinary  water  at  temperatures  between  10  and  20°C  may  be  ob¬ 
tained  from  standard  sources.7  The  thermal  expansivity  a  thus  obtained  is  equal  to 
0.00015°C'1 .  This  value  corresponds  to  a  Brunt-Vaisala  frequency  equal  to  0.0012  rad/s. 

In  the  GORWAK  right-handed  coordinate  system,  positive  x  is  the  direction  parallel 
to  the  sea  surface  away  from  the  stem.  Positive  y  is  starboard  of  a  person  standing  on  the 
deck  and  facing  the  bow  of  the  ship.  As  indicated  above,  the  positive  z  axis  points  to  the 
local  zenith.  In  order  to  simulate  the  measurements  of  Garrett  and  Smith2,  a  temperature 
decrease  of  1°C  for  every  10  m  of  increased  water  depth  is  assumed  throughout  the 
simulation.  As  the  dtpth  increases,  the  increments  of  both  T  and  z  are  negative.  Their 
ratio  is  positive  and  may  be  expressed  as 
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RESULTS 
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Making  use  of  a  bell-shaped  velocity  profile,  the  code  was  test-run  for  several  depths 
of  the  wake  axis  below  the  free  surface.  The  axis  of  the  wake  was  taken  from  a  depth  of 
20  m  to  about  10  m  below  the  surface.  Preliminary  results  showed  that,  for  wake  axis 
depths  shallower  than  10  in,  water  warmer  than  the  wannest  found  in  the  whole  profile 
appeared  right  at  the  surface,  clearly  a  meaningless  result. 

It  was  then  realized  that  GORWAK  performs  a  smoothing  (porosity)  interpolation 
both  at  the  water  surface  and  at  the  ocean  floor,  possibly  causing  some  anomalous  results 
at  these  boundaries.  These  suspicions  were  fully  confirmed.  By-pass  of  the  interpolation 
immediately  corrected  the  problem,  yielding  physically  meaningful  results  for  wake  depths 
as  shallow  as  2  m.  With  this  first  difficulty  successfully  overcome,  a  series  of  runs  was 
made  to  confirm  the  code  could  accurately  describe  the  behavior  of  the  macrowake. 

A  word  about  nomenclature:  we  noticed  that  the  data  we  received  were  labelled  dif¬ 
ferently  from  those  of  Garrett  and  Smith.2  In  an  effort  to  maintain  consistency  with  the 
published  literature,  we  inquired  as  to  the  discrepancy  and  were  informed  by  Smith  in  a 
private  communication  that  each  of  the  two  authors  has  adopted  his  own  system  of  wake 
identification.  The  one  used  in  their  paper  is  due  to  W.  D.  Garrett.  The  identification  of 
the  wakes  adopted  in  the  present  paper  is  the  one  used  by  P.  M.  Smith  and  found  on  the 
data  tapes  we  received  from  him. 

Program  input  data  are  listed  in  two  tables.  Table  1  contains  the  density  profile,  in 
oceanographic  sigma  units,  as  a  function  of  water  depth.  Table  2  contains  most  of  the 
other  important  parameters  used  in  the  same  simulation.  The  results  of  the  simulation  are 
extracted  from  computer  file  BL375.R6. 

The  data  in  Figs.  3  through  13  correspond  to  a  time  span  of  13.50  minutes.  They 
represent  the  time  evolution  within  a  vertical  plane  which  is  perpendicular  to  the  axis  of 
the  wake.  Swirl  is  given  a  value  of  zero.  The  problem  is  symmetrical  and  only  the  right 
half  of  the  plane  is  shown.  A  number  of  additional  results,  such  as  the  vorticity  and 
stream  functions,  are  listed  in  the  output  but  are  not  shown  here. 

The  GORWAK  output  is  expressed  in  terms  of  density  increments  which  must  be 
translated  into  temperature  changes.  The  density  of  sea-water  is  generally  a  function  of 
temperature,  pressure,  and  salinity.  The  persistence  of  the  macrowake  is  of  the  order  of  a 
few  hours.  The  time  constant  associated  with  temperature  changes  at  sea  is  of  the  order  of 
hours.  On  the  other  hand,  salinity  does  not  change  appreciably  for  days.  Hence,  for  the 
practical  purpose  of  this  study,  the  salinity  is  assumed  to  be  a  constant  equal  to  35  parts 
per  thousand.  This  number  is  a  realistic  average  value  of  ocean  salinity. s  Equation  4  can 
then  be  written  as  follows 


dg  =  (o  x  q)  dT.  (7) 

When  the  values  of  the  expansivity  and  the  density  near  the  surface  are  known,  each 
incremental  density  change  corresponds  to  one  value  of  the  temperature  difference.  Thus, 
a  temperature  difference  of  1  °C  corresponds  to  a  density  change  of  0,000172  or  1.72  x 
10'4.  Figures  3  through  13  are  drawn  to  represent  temperature  differences.  Note  that  an  in¬ 
frared  scanner  sees  only  the  vicinity  of  the  top  horizontal  line  in  the  plots,  which 
represents  the  free  surface. 
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Table  1.  Density  profile  for  macrowake  run  BL375.R6. 


Depth  (tn) 

0 

100 

J 

190 

200 

Density 

24.72 

26.44 

27.00 

27.10 

Note;  Density  values  are  in  oceanographic  urtits  of  sigma. 


Table  2.  Main  input  parameters  of  macrowake  run  BL375.R6. 


Variable  Name 

Description 

Value 

ACCUR 

Accuracy  parameter 

0.5 

BOTTOM 

Bottom  depth 

10000  in 

IDEPB  =  IDEPT 

Wake  depth 

3.7  m 

IPROFL 

Vertical  density  profile 

(measured  =  3) 

PARMIX 

Mixing  parameter 

0.8 

PORZ 

Smoothing  at  surface  and  bottom 

0 

SWIRL 

Amount  of  swirl 

0 

TOTTIM 

Converts  to  time  through  formula 

5 

WKRDFF 

Wake  radius 

10  m 

Following  the  P.  M.  Smith  wake  identification  scheme  already  mentioned,  we  in¬ 
vestigate  the  thermal  characteristics  of  Wake  7.  We  divide  the  contrast  field  into  nine  flat 
rectangular  elements.  The  length  of  each  rectangle  is  equal  to  the  distance  between  the 
profiles.9  The  width  of  the  rectangular  pieces  is  equal  to  the  width  of  the  corresponding 
profiles.  The  temperature  of  each  rectangular  element  is  the  average  of  the  corresponding 
profile. 

We  now  make  use  of  another  computer  program  called  SIRS.10  This  code  is  designed 
to  generate  the  thermal  characteristics  of  ships.  Vessels  under  study  are  divided  into  a 
number  of  flat  pieces.  In  order  to  model  the  wake,  we  substitute  the  rectangular  elements 
described  above  for  the  flat  elements  of  a  ship.  It  must  be  stressed  that  the  only  parts  of 
SIRS  that  are  used  in  this  application  are  the  atmospheric  propagation  and  the  detector 
modules.  The  temperature  of  the  fiat  pieces  is  computed  by  GORWAK  rather  than  SIRS. 
Using  parts  of  SIRS  as  we  did  was  dictated  purely  by  convenience.  Other  techniques  could 
have  been  used  and  might  have  performed  equally  well.  The  measured  and  theoretical 
£j  wake  results  are  compared  in  Table  3. 
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Table  3.  The  signal  from  Wake  7. 


MEASUREMENT 

THEORY 

Vert. 

nmi 

Range 

km 

Contrast  Radiant 
Intensity 

W/sr 

Contrast 

Irradiance 

W/cm2 

Contrast  Radiant 
Intensity 

W/sr 

Contrast 

Irradiance 

W/cm2 

1 

2 

— 0.298E+06 

-0.867E-05 

-0.225E  +  06 

-  0.655 E- 05 

5 

9 

-0.273E+06 

-0.318E-06 

-0.206E  F  06 

— 0.240E-06 

10 

19 

-0.270E+06 

-0.786E  '07 

— 0.204E  +  06 

— 0.594E-07 

20 

37 

-0.267E+06 

-0.194E-07 

— 0.202E  +  06 

-0.147E-07 

50 

93 

-0.267E+06 

-0.31  IE-08 

-0.202E  +  06 

-0.235E-08 

100 

185 

— 0.267E+06 

-0.778E-09 

— 0.202E+06 

— 0.588E-09 

200 

371 

-0.267E  +  06 

-0. 194E-09 

-0.202E  +  06 

—0. 147E-09 

300 

556 

-0.267E+06 

-0.864E-  10 

— 0.202E+06 

— 0.653E—  10 

SENSITIVITY  ANALYSIS 

A  sensitivity  analysis  of  the  GORWAK  computer  code  shows  how  strongly  each  of 
the  principal  factors  affects  the  results  and  which  particular  characteristics  of  the 
macrowake  are  affected.  Table  4  lists,  starting  from  the  left-hand  column,  the  physical 
quantities,  their  corresponding  variable  name,  the  program  default  value,  the  value  used  in 
the  calculation,  and  the  name  of  the  file  that  contains  the  computation  results.  Referring 
back  to  Table  2,  it  may  be  remembered  that  file  BL375.R6  contains  the  results  of  the 
closest  fit  simulation. 

The  variables  that  are  being  examined  arc:  the  macrowake  linear  velocity  (VWAKE), 
the  radial  steepness  (ANURAD),  the  radius  (WKRDFF),  and  the  distance  between  the  pro¬ 
peller  axis  and  the  free  water  surface  (IDEPB  =  1DEPT),  For  each  parameter,  the 
macrowake  profile  is  shown  at  three  distances  from  the  stem:  516  m,  1485  m,  and  2295 
m,  respectively.  The  VWAKE  results  are  contained  in  Figs.  14,  15,  and  16.  The 
ANURAD  variations  are  plotted  in  Figs.  17,  18,  and  19.  The  effect  of  varying  the 
macrowake  radius  WKRDFF  is  shown  in  Figs.  20,  21,  and  22.  Finally,  the  influence  of 
the  propeller  depth  is  displayed  in  Figs.  23,  24,  and  25.  In  Figs.  14  through  25,  line  plots 
indicate  measured  data.  They  are  shown  superimposed  on  best  fit  results. 

Another  program  variable  describes  swirl  (and  goes  by  the  name  SWIRL).  A  number 
of  simulation  runs  have  been  perfonned  to  examine  its  effects,  which  in  general  are  quite 
important,  on  the  macrowake.  However,  a  study  of  the  measured  data  does  not  show  any 
unsymmetrical  features  that  may  be  attributed  to  swirl.  For  this  reason,  as  Table  2  shows, 
the  SWIRL  variable  has  been  set  equal  to  zero  for  the  present  simulation. 
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Table  4.  Sensitivity  analysis  outline. 


Factor 

Wake  Linear 
Velocity 

Name 

VWAKE 

Wake  Radial 
Steepness 

ANURAD 

Wake  Radius 

! 

WKRDFF 

Propeller 

Depth 

1DEPB  -  IDEPT 

i 

1 _ 

Value  this  Run 


2  m 
3.7  m 
10  m 
20  m 


File  Name 


DM05.  R6 
BL375.R6 
Pi.374.R6 


BL375.R6 
DM04.  R6 
DM03.  R6 


BL377.R6 
BL375.R6 
BL381.R6 
DM08.  R6 


BLP2.R6 
BL375.R6 
DM10.R6 
DM12. R6 


SENSITIVITY  TO  THE  LINEAR  SPEED  OF  THE  MACROWAKE  (VWAKE) 

The  shape  of  the  macrowake  is  highly  sensitive  to  the  linear  speed.  Figures  15  (c) 
and  16  (e)  show  that  a  speed  of  10  em/s  Cuuscs  two  thermal  spikes  to  appear  at  a 
temperature  that  is  higher  than  any  other  found  in  the  thermal  profile.  This  nonphysical 
situation  is  a  limiting  ease.  A  careful  variation  of  VWAKE  reveals  that  a  speed  of  2  cm/s 
results  in  the  best  fit  with  the  measured  data. 

SENSITIVITY  TO  THE  SHAPE  OF  THE  VELOCITY  FIELD  (ANURAD  AND 
WKRDFF) 

The  ANURAD  variable  describes  how  sharply  the  velocity  falls  off  in  the  radial 
direction,  while  WKRDFF  is  an  effective  macrowake  radius,  measured  in  meters.  These 
two  variables  describe  the  shape  of  the  velocity  field.  Figures  19  and  22  show  that  tn  the 
far  field,  the  macrowake  profile  is  sensitive  to  changes  in  the  shape  of  the  velocity 
distribution. 

Starting  fror.:  an  examination  of  the  macrowake,  one  may  ponder  whether  it  is  possible 
to  draw  any  inferences  regarding  the  characteristics  of  the  ship  that  is  producing  it.  In  its 
present  shape,  the  GORWAK  c<xlc  stalls  from  a  field  velocity  distribution  input.  In  order  to 
find  an  answer  to  the  important  practical  problem  just  posed,  some  effort  will  be  required  if 
we  are  to  correlate  our  maciowake  results  with  the  design  parameters  of  the  ship. 
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SENSITIVITY  TO  THE  PROPELLER  DEPTH  (IDEPB) 

The  effect  on  the  macrovake  of  changing  the  propeller  depth  can  be  seen  in  Figs. 

23  through  25,  as  IDEPB  is  varied  from  2  m  to  20  m  beneath  the  free  water  surface.  The 
shape  of  the  maerowake  profile  is  quite  sensitive  to  this  variable.  During  the  simulation, 
IDEPB  is  set  to  a  value  of  3.7  m,  which  is  the  design  depth  of  the  Hayes  propellers  below 
the  sea  surface. 

Figure  23  depicts  a  macrowake  profile  caused  by  a  propeller  at  a  depth  of  20  m.  A 
small  scar  is  detectable  at  a  distance  of  516  m  behind  the  stem.  No  contrast  field  appears 
either  at  1485  m  (see  Fig.  24)  or  at  2295  m  (see  Fig.  25). 

DISCUSSION 

MODEL  ACCURACY 

The  overall  deviation  between  calculated  and  measured  values  in  Table  3  is  about 
24%.  Two  factors  conspire  to  make  the  agreement  somewhat  worse  than  it  ought  to  be. 

First,  the  vertical  thermal  profile  for  Wake  7  was  measured  at  least  two  days  before 
the  infrared  imagery  was  gathered.  Smith  says*  that  the  change  in  the  profile  over  these  48 
hours  cannot  be  estimated. 

Second,  the  GORWAK  program  explicitly  expects  the  local  thermal  background  to  be 
essentially  static  during  the  measurement,  The  underlying  assumption  is  that  two  spatial 
dimensions  plus  the  time  coordinate  are  equivalent  to  three  spatial  dimensions.  The  time 
variable  is  related  to  the  distance  behind  the  ship  via  a  constant  'inear  speed,  equal  to  4.71 
m/s  for  the  Hayes.  Several  clues  indicate  that  the  infrared  background  of  the  data  used  is 
active.  The  thermal  profiles  of  Fig.  26  are  distorted  to  the  right,  starting  some  2295  m 
behind  the  stem,  possibly  due  to  a  front  pushing  the  macrowake  to  the  right.  The 
photograph  representing  the  full  extent  of  the  measured  macrowake  clearly  shows  the 
presence  of  a  sharp  thermal  front  at  the  end  of  the  wake.11 

This  disturbance  notwithstanding,  the  agreement  between  measurement  and  theory  is 
still  acceptable.  Hence,  it  seems  that  the  program  can  be  relied  upon  to  yield  accurate 
results  even  in  the  presence  of  a  certain  level  of  disturbance,  as  leng  as  the  macrowake  is 
not  distorted  in  some  catastropnic  manner. 

THE  MACROWAKE  SHAPE  AND  SPATIAL  WIDTH 

Figure  26  compares  computational  data  with  measured  ones.  Plots  (a)  through  (f) 
have  only  one  minimum.  Plot  (g)  reveals  the  sudden  appearance  of  a  double  dip.  The 
computed  macrowake,  shown  on  the  same  graph,  goes  through  a  similar  double  minimum 
at  exactly  the  same  moment. 

The  computer-generated  spatial  width  of  the  macrowake,  drawn  in  Fig.  27  as  a 
continuous  line,  is  compared  with  measured  data  represented  by  black  dots.  Good  agree¬ 
ment  is  obtained  everywhere  except  for  the  last  point.  This  anomaly  may  be  explained  by 
the  sudden  appearance  of  the  thermal  front  already  mentioned. 


*  P.  M.  Smith,  in  a  private  communication. 
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GORWAK  simulation  results  have  also  been  compared  with  macrowake  width 
calculations  contained  in  a  study  of  the  interaction  between  microwave  radar  and  the  tur¬ 
bulent  wake.12  Some  of  the  computations  contained  in  that  report  are  applicable  to  the 
wake  of  a  small  (Class  1)  ship  of  5000  tons.  The  Hayes  is  in  the  same  weight  category, 
since  its  full  load  displacement  is  equal  to  3860  tons.1*  Swanson's  results  are  in  agreement 
with  the  values  displayed  in  Fig.  27. 

THE  MACROWAKE  GAP 

Next  is  a  test  of  the  near-field  behavior  of  the  theory.  Careful  examination  of  the  infrared 
imagery  reveals  that  in  some  instances  there  is  a  gap  at  the  onset  of  the  macrowake,  just 
behind  the  stem.  Figure  28  is  an  enhanced  photograph  of  the  surveying  snip  with  its  wake.  The 
gain  settings  on  the  infrared  scanner  made  the  signal  from  the  macrowake  just  discernible:  it  was 
estimated*  at  about  0.5°C.  The  length  of  the  gap,  measured  on  this  photograph,  is  275.5  m. 

Figure  29  is  a  black-and-white  photograph  of  a  false-color  coded  image  of  the  Hayes 
and  Wake  8.  The  dark  temperature  contour  closest  to  the  stem  corresponds  to  a 
temperature  difference  of  0.1  °C.  The  next  contour,  light  in  color,  represents  a  temperature 
difference  of  0.5°C.  For  comparison  with  the  previous  photograph,  the  contour  at  0.5°C  is 
considered.  The  distance  between  the  stem  and  the  0.5°C  contour,  measured  on  the 
photograph,  is  equal  to  269.1  m. 

The  calculated  maximum  temperature  contrast  of  the  macrowake  is  plotted  in  Fig.  30 
versus  the  distance  (the  contrast  is  negative).  A  contrast  of  0.5°C  corresponds  to  a 
distance  of  283  m  behind  the  ship.  The  measured  gap  length  of  the  Kane  deviates  from  the 
theory  by  2.7%,  the  one  on  the  Hayes  differs  by  4.9%. 

CONCLUSIONS 

Using  a  generic  bell-shaped  velocity  distribution  and  no  swirl,  we  have  been  able  to 
simulate  the  infrared  signal  generated  by  the  macrowake  of  a  twin-hull  surface  ship,  the 
Hayes ,  running  only  on  its  starboard  engine.  The  contrast  radiant  intensity  of  the 
macrowake,  the  evolution  and  extent  of  the  spatial  width,  and  the  gap  length  all  Fit  the 
measurements.  Impressive  as  the  results  are,  we  do  not  yet  feel  justified  in  making  sweep¬ 
ing  statements  that  apply  to  all  surface  ship  wakes.  At  present,  it  must  be  remembered  that 
only  one  macrowake  measurement  has  been  analyzed.  The  preceding  results  raise  an  in¬ 
teresting  question.  Figure  13  shows  the  orescnce  of  large  finger-like  components  of  the 
wake  inside  the  water.  Does  such  a  flow  pattern  really  occur  in  the  wake?  The  only 
evidence  wc  now  have  of  that  possibility  appears  in  Reference  14.  It  is  clear  that  careful 
detailed  measurements  are  needed  inside  the  water,  down  to  depths  of  several  tens  of 
meters. 

The  experience  gained  in  this  project  leads  to  the  following  criteria  as  a  basis  for 
making  good  quality  infrared  macrowake  measurements.  To  be  useful,  infrared  imagery 
should  preferably 

1.  Be  taken  in  a  down-looking  geometry,  from  an  airborne  platform  flying  along  the  axis 
of  the  wake  and  using  a  top  quality  calibrated  stable  infrared  viewer. 

2.  Be  accompanied  by  a  vertical  air  and  water  temperature  profile  measured  as  closely  as 
possible  in  time  and  space  with  the  infrared  data. 


*  P.  M.  Smith,  in  a  private  communication. 
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3.  Include  thermal  profiles  measured  across  the  macrowake  at  several  distances  behind  the 
ship  and  at  equal  depths  up  to  about  50  m  beneath  the  sea  surface. 

4.  Include  other  data  as  well,  such  as  information  about  the  ship  producing  the  wake, 
other  accompanying  wakes,  ship  hull  design  specifications,  measurements  of  the  flow 
field  around  the  hull,  propeller  design,  near-field  propeller  flow  in  front  of  and  behind 
the  propeller,  ship  engine  power  specifications,  cooling  and  discharge  temperatures, 
past  thermal  ship  history  and  movements,  Oceanic  currents,  and  possible  presence  of 
older  wakes  in  the  area. 

Work  on  the  other  available  infrared  imagery  has  been  started.  It  is  unfortunate  that 
we  lack  a  measurement  of  the  exact  velocity  distribution  associated  with  the  macrowakes 
mentioned  in  the  present  analysis.  We  also  lack  careful  macrowake  measurements  on  other 
ship  types  and  classes,  be  they  full-scale  or  to-scale  models.  We  hope  the  present  effort 
will  encourage  more  research  in  this  area  and  emphasize  the  need  for  careful  gathering  of 
some  urgently  needed  measurements. 
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The  infrared  surface  ship  macrowake  problem.  (The  measurement 
instrumentation  is  flown  on  board  the  aircraft.) 


Fig.  2.  Block  diagram  of  the  theoretical  scheme. 
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Fig.  5.  The  calculated  tem|»erature  difference  contours  at  2.70  min.  (Contours  are 
displayed  in  a  vertical  section  through  the  macrowake.  The  pattern  is 
symmetrical  about  the  left-hand  side  vertical  axis.) 
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Fig.  6.  The  calculated  temperature  difference  contours  at  4.08  min.  (Contours  are 
displayed  in  a  vertical  section  through  the  macrowake.  The  pattern  is 
symmetrical  about  the  left-hand  side  vertical  axis.) 
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Fig.  II.  The  calculated  temperature  difference  contours  at  10.80  inin.  (Contours  arc 
displayed  in  a  vertical  section  through  the  macrowake.  The  pattern  is 
symmetrical  about  the  left-hand  side  vertical  axis.) 
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Fig.  12.  The  calculated  temperature  difference  contours  at  12.18  min.  (Contours  are 
displayed  in  a  vertical  section  through  the  macrowake.  The  pattern  is 
symmetrical  about  the  left-hand  side  vertical  axis.) 
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Fig.  13.  The  calculated  temperature  difference  contours  at  13.50  min.  (Contours  are 
displayed  in  a  vertical  section  through  the  macrowake.  The  pattern  is 
symmetrical  about  the  left-hand  side  vertical  axis.) 
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Fin.  14.  The  macrowakc  profile  at  516  m  with  VWAKE  equal  to:  (a)  0  cm/s, 
(h)  2  cm/s,  (c)  10  cm/s. 
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Fig.  15.  The  macrowake  profile  at  1485  in  with  VVVAKE  equal  to:  (a)  0  cm/s, 
(b)  2  cm/s,  (c)  10  cm/s. 
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Fig.  19.  The  macrowake  profile  at  2295  m  with  ANURAD  equal  to:  (a)  2,  (b)  4, 
(c)  10. 
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Fij*.  20.  The  macrowake  profile  at  516  m  with  WKRDFF  equal  to:  (a)  5  m,  (b)  10  m 
(c)  15  m,  (d)  20  m. 
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(d)  NO  MACRO-WAKE 

0  10  20  30  40  50  60  7  0  80  90  100 
M 

The  macrowake  profile  at  2.295  ni  and  propeller  depth  of:  (a)  2  m 
(b)  3.7  ni.  (c)  10  ni,  (d)  20  m  (no  macrowake  appears). 
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Fig.  26.  Comparison  of  the  Wake  7  measured  profiles  (drawn  with  a  continuous  line) 
with  theoretical  calculations  (drawn  with  a  hash  mark).  The  corresponding 
times  or  distances  from  the  stern  of  the  Hayes  are:  (a)  1.1  min.  or  311  m 
(b)  1.83  min.  or  516  m  (c)  2.47  min.  or  697  m  (d)  3.93  min.  or  1110  m 
(e)  5.25  min.  or  1485  m  (D  6.69  min.  or  1890  m  (g)  8.12  min.  or  2295  in 
(h)  9.54  min.  or  2695  m  (i)  10,86  min.  or  3070  m. 
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Fig.  27.  The  spatial  width  measured  across  Wake  7  on  the  free  surface  (drawn  with  a 
continuous  line),  compared  with  the  calculated  values  (shown  as  full  dots). 
The  macrowake  was  generated  by  the  Hayes. 
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Fig,  28.  A  digitally  enhanced  photograph  shows  the  Kane  and  part  of  the  macrowake 
it  generates.  (Notice  the  gap  in  this  8  to  12  micrometer  band  image.) 
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Fig.  30.  Plot  of  the  calculated  near-fleld  temperature  contrast  of  the  macrowake,  as  a 
function  of  the  distance  from  the  stern.  (Numerical  values  are  taken  from 
file  BL375.R6.) 
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